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Abstract. We considered the lattice electroweak theory at realistic values of a and 
9w and for large values of the Higgs mass. We investigated numerically the properties 
of topological objects that are identified with quantum Nambu monopoles. We have 
found that the action density near the Nambu monopole worldlines exceeds the density 
averaged over the lattice in the physical region of the phase diagram. Moreover, their 
percolation probability is found to be an order parameter for the transition between 
the symmetric and the broken phases. Therefore, these monopoles indeed appear 
as real physical objects. However, we have found that their density on the lattice 
increases with increasing ultraviolet cutoff. Thus we conclude, that the conventional 
lattice electroweak theory is not able to predict the density of Nambu monopoles. This 
means that the description of Nambu monopole physics based on the lattice Weinberg 
- Salam model with finite ultraviolet cutoff is incomplete. We expect that the correct 
description may be obtained only within the lattice theory that involves the description 
of TeV - scale physics. 



PACS numbers: 12.15.-y, 11.15.Ha, 12.10.Dm 
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1. Introduction 

The electroweak theory does not contain topologically stable monopole-like objects. 
However, certain unstable objects of topological nature still exist in this theory. One 
of the examples is the so-called Nambu monopole[^. It must be connected by the so- 
called Z string with the corresponding antimonopole. The Z string has nonzero tension. 
Therefore only the monopole - antimonopole bound state may appear as an observable 
object. The mass of the Nambu monopole (realized as a classical field configuration) 
was estimated to be of the order of several TeV. This is not far from the energies that 
may be achieved by modern colliders, in particular, the LHC Thus one may suppose, 
that an indication of its existence may be detected in the near future. 

Nambu monopoles are not described by means of a perturbation expansion around 
the trivial vacuum background. Therefore, nonperturbative methods should be used in 
order to investigate their physics. Lattice methods seem to be one of the ways to deal 
with Nambu monopoles. It should be stressed that the mass of the Nambu monopole is 
close to the energy scale, where (as commonly believed) the Standard Model does not 
work[2],[3]. This creates an additional difficulty while considering the problem. 

Qualitative lattice investigations of the properties of Nambu monopoles in the 
Standard Model have been performed both at zero and finite temperature in the 
unphysical region of large coupling constants [HElEllTj. Nambu monopoles were found 
to be condensed in the symmetric phase of lattice theory (and above the electroweak 
transition in the finite temperature theory). In the present paper we continue this 
investigation for realistic values of the renormalized coupling constants {a ~ 1/128 and 
6w = 7r/6) within the zero temperature theory. It should be stressed that originally 
Nambu monopoles were defined as classical objects^. Therefore there could appear 
several lattice definitions of quantum Nambu monopoles. In this paper we discuss two of 
them and investigate the difference between the positions of the corresponding monopole 
trajectories. 

The numerical investigation of Gauge-Higgs models has a long history. First studies 
of nonabelian Gauge-Higgs systems were performed in eighties and were devoted mainly 
to the investigation of SU (2) Gauge-Higgs model at zero temperature (see, for example, 
P, El [10, [m [12] and references therein). In these studies the general phase structure of 
the model was established. At small values of the scalar self coupling A the first order 
phase transition between the physical Higgs phase and the unphysical symmetric phase 
of the model was found. The phase transition becomes weaker as A is increased. Using 
the weak coupling expansion [1.3| it was found that the maximal value of the ultraviolet 
cutoff in the theory is achieved at infinite A. At A cxo also the absolute upper bound 
on the Higgs boson mass is achieved which was found to be of the order of lOM^ (see, 
for example, [IDl [HI [12] )• An early study of the SU{2) ® U{1) Gauge-Higgs model was 
performed in [H], where the phase diagram of the model was given. In general, the 
description of the phase diagram is similar to that of given in the present paper (see 
Fig. [H and discussion in section 5.1). 
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Table 1. The values of the cutoff used in some selected lattice studies of the SU{2) 
Gauge - Higgs model. 



Reference 


ultraviolet Cutoff (GeV) 


Mh (GeV) 




[15J 


140 (space direction) 570 (time direction) 


80 




16J 


280 (time direction) 


80 


[17J 


280 


34 




[18J 


110 


16 


il9j 


90 (space direction) 350 (time direction) 


34 


[20] 


280 


48 


[21] 


140 


35 


[22J 


280 


20 , 50 


[23J 


190 


50 




24J 


260 


57-85 


[25J 


200 - 300 


47 - 108 


[lOj 


400 


480 




[llj 


330 - 470 


280 - 720 (both A = oo and finite A) 


[12J 


250 - 470 


720 (both A = oo and finite A) 



The next step in the numerical investigation of Gauge- Higgs models was motivated 
by an attempt to explore the physics of finite temperature electroweak phase transition 
and its relation to cosmology, in particular, to the problem of baryon asymmetry (see, 
for example, [l5l[l6l[I3[l8l[19l[20l[2ll[2a[23l[H[25l references therein) . One of 

the achievements of these studies was that the endpoint of the phase transition line (in 
the T — Mh plane) was found. It was found that at the experimentally allowed values 
of the Higgs mass the transition in the standard electroweak theory is a crossover. 

In AD lattice studies of SU{2) Gauge-Higgs theory the scale is fixed by the value of 
the W^-boson mass. Namely, the mass measured in lattice units is My/ = a x 80 Gev, 
where a is the value of the lattice spacing. The ultraviolet cutoff is A = ^. In Tabled] 
we summarize the data on the ultraviolet cutoff used in selected lattice studies of the 
SU (2) Gauge-Higgs model. In all these studies the coupling constant corresponding to 
the fine structure constant of the Weinberg-Salam model is around a ~ 1/100. The 
potential for the Higgs boson was used with different values of A. 

The numerical investigation of the monopole-like [27] and string - like topological 
objects in Abelian Gauge-Higgs systems has also a long history (see, for example, |28j). 
For an early study of the percolation of monopole currents see [29j. It has been found 
that the monopole density drops sharply to zero in the Coulomb phase of the model 
while in the confining phase it is nonzero. The vortex density is nonzero in the Coulomb 
phase and drops to zero in the Higgs phase. The situation with topological defects in 
nonabelian Gauge-Higgs systems is in general similar to that of the Abelian ones. The 
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first numerical investigation of topological defects in 3D lattice SU{2) Gauge-Higgs 
model (corresponding to finite temperatures) was performed in [30]. It was found that 
the percolation of Z-vortices is the order parameter for the electroweak transition. The 
conjecture that Nambu-monopole percolation feels the electroweak transition was made 
first in |31j . 

In our earlier papers we considered the appearance of an additional discrete 
symmetry in the fermion sector of the Standard Model[H 13 [6], [7]. This additional 
symmetry allows to define the Standard Model with the gauge group SU (3) x SU (2) x 
U{l)/Z, where Z is equal to Zq, or to one of its subgroups: or Z2. It is worth 
mentioning that it has been recognized much earlier that the Standard Model appears 
with the gauge group SU{3) x SU (2) x f/(l)/ Zq as a result of the spontaneous breakdown 
in the SU{5) unified model|32j. Independently, the Zq symmetry in the Higgs sector of 
the Standard Model was considered in [33] . 

In the present paper we use two lattice realizations of the electroweak theory: 
with the gauge groups SU{2) x U{l)/Z2, and SU{2) x f/(l), respectively. The 
SU (2) xU{l)/Z2 model should be the part of the Standard Model with the gauge group 
SU{3) X SU{2) X U{l)/Z2 or SU{3) x SU{2) x U{1)/Zg, while the SU{2) x f/(l) model 
could be the part of the Standard Model with the gauge group SU{3) x SU (2) x f/(l)/Z3 
or SU{3) X SU{2) x f/(l). We comment on the difference between the two lattice models, 
which appears in our numerical research in the unphysical region with a large coupling 
constant. 

The paper is organized as follows. In Sect. 2 we describe the lattice models under 
consideration. In Sect. 3 we present the definition and the main properties of quantum 
Nambu monopoles. Sect. 4 contains our description of the quantities to be measured. In 
Sect. 5 we report our main numerical results, while in Sect. 6 we discuss the difference 
between the two versions of the lattice electroweak model. The final section contains 
our conclusions. 

2. Lattice models under investigation 

In this section we describe the lattice models under consideration. We do not consider 
the color sector of the Standard Model. Therefore, we are left with two possibilities: the 
gauge groups SU (2) x U {1)/Z2, and SU (2) x f/(l). We also neglect dynamical fermions. 
In both cases we use the following lattice variables: 
1. The gauge field U = (U, 9), where 




(1) 



are realized as link variables. 
2. A scalar doublet 



a = 1,2. 



(2) 
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The potential for the scalar field is considered in its simplest form [S| in the London 
limit, i.e., in the limit of infinite bare Higgs mass. In the lattice study this does not 
mean, however, that the physical Higgs mass is infinitepj. Instead we expect only that 
it should not be less than the inverse lattice spacing. This is indeed confirmed via direct 
calculation. From the very beginning we fix the unitary gauge $i = y/j, $2 = 0. 

For the case of the SU{2) x U{l)/Z2 symmetric model we chose the action of the 

form 

Sg = P {{I - I T^r Up COS 9p) + 1(1- COS 29p)) + 

plaquettes 

+ij:i^-Muiy^^)), (3) 

xy 

where the plaquette variables are defined as Up = U^yUyzU^^U*^, and 9p = 6.j.y + Oy^ — 
Owz — Gxw for the plaquette composed of the vertices z, w. 

For the case of the conventional SU{2) x U{1) symmetric model we use the action 

Sg = P J2 {{1 - Up) + 3(1 -COS Op)) + 

plaquettes 

+ lT.i^-Re{Uiy^y)). (4) 

xy 

In both cases the bare Weinberg angle is 6w = vr/6, which is close to its experimental 
value. The renormalized Weinberg angle is to be calculated through the ratio of the 
lattice masses: cos9w = Mw/Mz- The bare fine structure constant a is expressed 
through j3 as a = l/AnjS. However, the renormalized coupling extracted from the 
potential for infinitely heavy fermions differs from this simple expression, as will be 
shown in the next sections. The physical meaning of the constant 7 is that it is equal 
to the square of the vacuum value of $1. 

The following variables are considered as creating a Z boson and a W boson, 
respectively: 

Z^y =Z^ =sm[ArgUll+e,yl 

Wxy = W>i = Uiy-y. (5) 

Here, n represents the direction {xy). 

After fixing the unitary gauge the electromagnetic U{1) symmetry remains: 

Uxy ^ 9xUxy9y^ 

dxy dxy - 0Lyl2 + (6) 

where = diag(e*°^/^, e~*"^/^). 

In the unitary gauge there is also a U{1) lattice gauge field, which is defined as 

Axy = = [-Argf/i; + e,y] mod 277, (7) 

The fields A, Z, and W transform as follows: 



7 N 7 
^xy ' ^xyi 

Wxy^Wxye-'""^. (8) 
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It should be mentioned that the field A cannot be treated as the usual electromagnetic 
field, because the set of variables A, Z, and W do not diagonalize the kinetic part of the 
pure gauge action in its naive continuum limit. In our lattice model the electromagnetic 
field Aem should be defined as 

Aem = A + Z' -2 sin^ ewZ', (9) 

where Z' = [AigU^l + 6'^j^]mod27r. 

3. Nambu monopoles 

First, we define the continuum electroweak fields as they appear in the Weinberg-Salam 
model in the way appropriate for the topological consideration. Namely, after fixing the 
unitary gauge $i = const., $2 = 0, where $ is the scalar field of the electroweak theory, 
the Z-boson field Z'^ and electromagnetic field are defined as 

A^M = 25^-2 sin^ OwZ^, (10) 

where and are the corresponding SU{2) and U{1) gauge fields of the Standard 
Model. 

Nambu monopoles are defined as the endpoints of the so-called Z-string [Ij. The 
Z-string is the classical field configuration that represents an unstable object, which 
is characterized by the magnetic flux extracted from the Z-boson field. Namely, for a 
small contour C winding around the Z - string one should have 

J Z^dzf" ~ 27r; j A'^y^dx^ ^ 0; J B^'dx^' ~ 27rsin2^iy (H) 

The string terminates at the position of the Nambu monopole. The hypercharge flux 
is supposed to be conserved at that point^. Therefore, a Nambu monopole carries 
electromagnetic flux 47rsin^ 9w- The size of Nambu monopoles was estimated [1] to be 
of the order of the inverse Z-boson mass, while its mass should be of the order of a few 
TeV. According to [1] Nambu monopoles may appear only in the form of a bound state 
of a monopole-antimonopole pair. 

In lattice theory the classical solution corresponding to a Z-string should be formed 
around the 2-dimensional topological defect which is represented by the integer- valued 
fleld deflned on the dual lattice 

J: = ^ {[dZ'U^2. - dZ') (12) 

(Here we used the notations of differential forms on the lattice. For a deflnition of those 
notations see, for example, [3lj.) Therefore, S can be treated as the worldsheet of a 
quantum Z-stringfSTl [TJ [30] . 

I On the classical level the monopole - like topological objects with nontrivial hypercharge flux have 
an infinite self energy. At the quantum level such objects are present in the lattice theory at the finite 
values of lattice spacing. However, in the physically interesting region of lattice coupling constant 
/3 ~ 15 the density of these objects vanishes (see Section 6 of the present paper). 
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Then, worldlines of quantum Nambu monopoles appear as the boundary of the 
Z-string worldsheet: 

Jz = 5S (13) 

It has been mentioned in the previous section that our lattice models become U{1) gauge 
models after fixing the unitary gauge. The corresponding compact U{1) gauge field is 
given by Eq. (JTj). Therefore one may try to extract monopole trajectories directly from 
A. Actually this was done in our earlier papers [H |5[ [6], [7j . The monopole current is 
given by 

JA = ^*di[dA]mod2Tr) (14) 

Both jz and Ja represent objects carrying magnetic charge. Therefore it would be 
instructive to reveal the correspondence between them. We have 

A = [-Z' + 2e]mod27r. (15) 

In continuum notation this would be 

= + 23", (16) 

where B is the hypercharge field. Its strength is divergenceless. As a result in continuum 
theory the net Z flux emanating from the center of the monopole is equal to the net A 
flux with the opposite sign. (Both A and Z are undefined inside the monopole.) This 
means that in the continuum limit the position of the Nambu monopole must coincide 
with the position of the antimonopole extracted from the field A. Therefore, one can 
consider Eq. (HM as another definition of a quantum Nambu monopole. It is interesting 
that the definition (fT^ is not directly related to any observable string, as the Dirac 
string connecting the corresponding lattice monopoles is invisible. 



4. Quantities to be measured 
4.I. Evaluation of the lattice spacing 

The physical scale is given in our lattice theory by the value of the Z-boson mass 
Ml*"^" ~ 90 GeV. Therefore the lattice spacing is evaluated to be a ~ [90 GeV] ^Mz, 
where Mz is the Z boson mass in lattice units. 

In order to evaluate the mass of the Z-boson we use the correlator [9]: 

Y.iY, Z^Z^) ~ e-*^^l^«-^«l + e-^-^^(^-l^'«-?'«l), (17) 

x,y /I 

Here the summation J2x,y is over the three "space" components of the four - vectors x 
and y while XQ,yQ denote their "time" components. L is the lattice length in the "time" 
direction. 

It is worth mentioning, that in the Z-boson channel many photon states also exist. 
The mass of the corresponding state on the finite lattice we used is, however, larger 
than that of the Z - boson itself. For example, on the lattice 16'^ x 24 the minimal mass 
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of the 3 - photon state is M^^ = 2^ + 1.5. Moreover, from the point of view of 
perturbation theory this state appears in the correlator (|T71) through a virtual loop and 
is suppressed by the factor a^. 

4-2. The Higgs boson mass. 

The Higgs boson mass in lattice units is measured using the correlator 

Y^iiH^Hy) - {Hf) ~ e-*^^!""-^'"! + e-^^^^^-l^^-^^oD, (18) 

where H is the Higgs boson creation operator. 

We used three different operators that create Higgs bosons: 

H^ = Y,\W,y\\ (19) 

y 

H. = T.Zly (20) 



xy 

y 



and 



H^ = Y.Re{Uiy^^) (21) 
y 

Here is defined at the site the sum ^y is over its neighboring sites y. 
4.3. The renormalized coupling 

The bare constant a = e'^/An (where e is the electric charge) can be easily calculated 
in our lattice model. It is found to be equal to l/(47r/3). Therefore, its physical value 
a{Mz) ~ 1/128 could be achieved at values of f3 in the vicinity of 10. This naive guess 
is, however, to be corrected by the calculation of the renormalized coupling constant or. 
We perform this calculation using the potential for infinitely heavy external fermions. 
We consider Wilson loops for the right-handed external leptons: 

Wi^p,(/) = (ReH(.,),,e2^^-). (22) 

Here I denotes a closed contour on the lattice. We consider the following quantity 
constructed from the rectangular Wilson loop of size r x t: 

V(r) = log lim ^(^ x^) (23) 
^ ^ ^ t^c»>V(r X (t + 1)) ^ ^ 

Owing to the exchange of virtual photons at large distances we expect the appearance 

of the Coulomb interaction 

V(r) = -— + const. (24) 
r 

It should be mentioned here, that in order to extract the renormalized value of a one 
may apply to V the fit obtained using the Coulomb interaction in momentum space. 
The lattice Fourier transform then gives 

V(r) = — a/jW(r) + const, 

Uir) = 4-.y ^ (25) 

' L^f^^ + sm2p2/2 + sm2p3/2 ^ ' 
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Here L is the lattice size, Pi = ^ki,ki = 0, L — 1. On large enough lattices at r << L 
both definitions approach each other. For example, for L = 75, r G [1, 10] the linear fit 
to the dependence U{r) on ^ gives U{r) ~ 0.97/r — 0.18 while for L = 100, r G [1, 10] 
the fit is U{r) ~ 0.997/r — 0.155. However, on lattices of the sizes we used the difference 
is important. Say, on the lattice 24^ the fit is U{r) ~ 0.82/r — 0.35 (for r G [1,5]). 
Thus, the values of the renormalized extracted from flM|) and fl25|) are significantly 
different from each other. Any of the two ways, fl24|) or fl25|) . may be considered as the 
definition of the renormalized a on the finite lattice. And there is no particular reason 
to prefer the potential defined using the lattice Fourier transform of the Coulomb law 
in momentum space. Actually, our study shows that the single 1/r fit approximates 
V much better. Therefore, we used it to extract a/j. This should be compared with 
the results of [12], where for similar reasons the single e~^^/r fit (instead of the lattice 
Yukawa fit) was used in order to determine the renormalized coupling constant in the 
SU (2) Gauge-Higgs model. However, the fact that both definitions give values that 
differ from each other shows the limitation on the interpretation of our results and the 
importance of finite volume effects in our research. 



4.4- Namhu monopole density and percolation probability 

According to Eqs. ( |T3l [Ml) the worldlines of the quantum Nambu monopoles could be 
extracted from the field configurations in two ways: 

jz = 5S = —*d{[dZ']mod2'K) (26) 
2tt 



and 



= 5S = ^*d{[dA]mod27i). (27) 
2tt 

The monopole density is defined as 

(X^links |jlink| \ fno\ 
/' ^^^^ 

where L is the lattice size (in lattice units). 

In order to investigate the condensation of monopoles we use the percolation 
probability H(A). It is the probability that two infinitely distant points are connected 
by a monopole cluster (for more details of the definition see, for example, |35j). 

Both —Ja and +jz describe the same physical object. However, this object may 
have a size that is larger than one lattice spacing. That's why the two different ways 
to extract the monopole trajectory may give different currents. The difference between 
the two currents is jz — (— Ja) = jz + Ja- Therefore, the density of Ja + jz measures 
the degree of how Ja differs from —jz- In order to investigate the difference between 
the two definitions of Nambu monopole currents we use the quantity p(j^ + jz), that is 
constructed using the current jz + jA as in ( l28l) . 
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4-5. Action density near monopole trajectories 

The monopole worldline lives on the dual lattice. Each point of the worldline is 
surrounded by a three - dimensional hypercube of the original lattice. We measure 
the plaquette part of the action S^°^ on the plaquettes that belong to those three- 
dimensional hypercubes (normalized by the number of such plaquettes). The excess 
of the plaquette action near monopole worldlines over the mean plaquette part of the 
action Sp is denoted by 

A5p=^(^r°-^p)- (29) 

Op 

Very roughly ASp can be considered as measuring the magnetic energy (both SU{2) 
and U{1)), which is carried by Nambu monopoles. 

We also measure Sf^"^, which is the part of the action 5*™°"^ on the links of the 
original lattice that connect vertices of the two incident three-dimensional hypercubes 
mentioned above. The excess of this link action near monopole worldlines over the mean 
link part of the action Si is denoted by 

ASi = Usr''-Si). (30) 

For the simplicity of the calculations we use only one of the 8 links that connect incident 
hypercubes. 

4-6. Hypercharge monopoles 

In addition to the Nambu monopoles we also investigated the behavior of some objects 
which are called hypercharge monopoles. Their worldlines are extracted from the 
hypercharge field 6 in the following way: 

jy = —*d{[d2e]mod2n) (31) 

271 

We also define their density according to the expression fl25]) . 

Actually, in the naive continuum limit hypercharge monopoles would have an 
infinite energy. They may appear only if one takes into account the finiteness of the 
ultraviolet cutoff. It occurs that they are mainly of interest in the strong coupling region, 
where the two considered lattice models appear to behave differently. 

5. Numerical results 
5.1. Phase diagram 

The phase diagrams of the two models under consideration are presented in Fig. [H At 
small values of the coupling constants the model with the gauge group SU (2) U{l)/Z2 
has already been investigated in our earlier paper [5]. The model with the gauge 
group SU{2) ® U{1) was investigated in the paper [H] (also at small values of the 
coupling constants). The dashed vertical line represents the phase transition in the 
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SU{2) ® t/(l)-symmetric model (we call it further Model A). This is the confinement- 
deconfinement phase transition corresponding to the U{1) constituents of the model. 
The same transition for the SU{2) ® f/(l)/2'2-symmetric model (we call it model B) is 
represented by the solid vertical line. The dot-dashed horizontal line corresponds to the 
transition between the broken and symmetric phases of model A. The solid horizontal 
line represents the same transition in model B. Interestingly, in the SU{2) ® U{l)/Z2 
model both transition lines meet, forming a triple point. Much attention was paid to 
this fact in [5]. 

So, in both models there are three phases. The first one is situated in the left- 
hand side of the phase diagram. In this phase there are confinement-like forces both 
between the right-handed and the left-handed external fermions. However, due to the 
presence of the charged scalar field the string connecting external fermions is broken 
and the confining forces disappear at a certain distance (see, for example, \5\). In this 
phase both Nambu monopoles and hypercharge monopoles are condensed. The second 
phase is situated below the horizontal phase transition line and right to the vertical 
phase transition line. In this phase the confining forces are observed only between 
the left-handed fermions. The hypercharge monopoles are not condensed in this phase 
and their density falls sharply. For the detailed description of different phases of the 
SU{2) ® f/(l)/Z2-symmetric model (including the properties of topological objects and 
static potentials) at small values of (3 see [5]. 

Real physics is commonly believed to be achieved within the phases of the two 
models situated in the right upper corner of Fig. 1. In this phase neither Nambu 
monopoles nor hypercharge monopoles are condensed. The confining-like forces are 
not observed here both between right-handed and left-handed fermions. The double- 
dotted-dashed vertical line on the right-hand side of the diagram represents the line, 
where the renormalized a is constant and equal to 1/128. In order to draw the phase- 
transition lines at small values of (3 we use the results of [5] and These data have 
been checked using the observables listed in the previous section. In particular, the 
density of hypercharge monopoles appears to be very sensitive to the U{1) confinement- 
deconfinement phase transition, while the density and percolation probability of Nambu 
monopoles feel the transition between the broken and symmetric phases. The position 
of the horizonal line for large values of (5 is obtained using mainly the percolation 
probability for the Nambu monopoles (which has been found to be an order parameter 
of the corresponding transition in the lattice Standard Model, with SU (3) constituents 
included [7]). This position corresponds also to the maximum of the susceptibility 
X = {H"^) — {H)"^. (See Fig. [6l We calculated the susceptibility using both definitions, 
Eqs. (IHand ([20]) for H.) 

All simulations were performed on lattices of sizes 8^ and 16^. Several points were 
checked using a lattice 24^. In general we found no significant difference between the 
mentioned lattice sizes. 
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5.2. Renormalized masses and couplings 

In the region (3 G (10,20), 7 G (1,2) we found no difference between the two versions 
of lattice electroweak theory. Therefore, we omit mentioning to what particular model 
the considered quantity belongs in this region of coupling constants. 

The Z-boson masses is found to change very slowly with the variation of The 
dependence on 7 seems to be stronger. Mz in lattice units grows with the increase of 7. 

For the calculation of the Z-boson mass we used lattices of sizes: 6^ x 12, 8^ x 16, 
12^ X 24, and 16'^ x 24. The dependence of the Z-boson correlator Eq. ( [T7|) on r = Xq — 2/0 
is presented in Fig. 2 for 7 = 1, /3 = 15 and a lattice 16^ x 24. From this plot we extract 
the mass of the Z-boson to be 0.22 ± 0.02. It is important to notice that we did not 
find any dependence of Mz on lattice size. 

At /? = 15 we localize the position of the transition between the symmetric and 
broken phases of the model at 7c = 0.92 ± 0.02. The measured Z - boson mass for 
7 G [0.9, 0.94] is 0.21 ± 0.02. So, we evaluate the Z - boson mass at the transition point 
to be Mz = 0.21 ± 0.02. The correspondent value of the ultraviolet cutoff (the inverse 
lattice spacing) is A = 90 GeV/Mz = 430 ± 40 GeV. 

As for the Higgs boson mass, due to the insufficient statistics we cannot extract 
Mh from our data with reasonable accuracy. According to our (very rough) estimate 
at /3 = 15, 7 = 1 we have Mh/Mz ~ 9 ± 2. This estimate is in agreement with the 
investigation of the SU (2) Higgs model [TOl |TT1 |T2] performed near the transition point 
for the London limit of the Higgs potential and realistic (3. Actually, as in [10] we made 
our estimate based on the consideration of the correlator for small "time" separation 
(< 3). 

It was found in [12] that at larger distances a second mass parameter close to 2Mw 
contributes to the correlator. In our study the accuracy of measurements does not 
allow us to extract information from the H-H correlator for "time" separations > 4. 
Therefore, we do not see the signal of the two gauge-boson bound state. However, we do 
not exclude that it would appear if more statistics is collected and "time" separations 
> 4 are explored. 

In [12] in order to evaluate the Higgs boson mass in this situation the value of 
mass ~ 2Mw was considered as the mass of the bound state of the two gauge bosons, 
and only the first mass in the given channel was interpreted as the Higgs boson mass. 
Actually, we do not see any reason to do so, and guess that the second mass in this 
channel may serve as the Higgs boson mass. However, this question must be investigated 
separately. Therefore we only claim here that the Higgs boson mass for our choice of 
initial parameters is larger than about 2Mw ^ 160 GeV, that could be the lowest mass 
in the given channel. 

The renormalized coupling constant a is found to be close to the realistic value 
a{Mz) = 1/128 along the line represented in Fig. 1. In Fig. 3 the dependence of the 
potential for infinitely heavy right-handed leptons on 1/r is shown for 7 = l,/5 = 15. 
The renormalized is extracted from this dependence. Actually, a linear dependence 
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is observed for r G [1,6]. Therefore we treat this constant as aR{r/a) ~ a/j(M^). In 
Fig. 4 we exhibit the dependence of l/a^ on (3 for fixed 7 = 1. Here it should be 
mentioned that according to the subsection IV. C of the present paper, the definition of 
the renormalized a on the lattices of finite size suffers from lattice artifacts. 

5.3. Namhu monopoles 

We used both definitions of Nambu monopoles given in section 3. In Fig. 5 we 
show their density and percolation probability as a function of 7 along the line 
of constant renormalized = 1/128. Interestingly, the density and percolation 
probability coincide here for the two mentioned definitions of Nambu monopole while 
the precise position of monopole trajectories differ by about 30%, i.e., we found that 
2p(jA + jz)/ {pUa) + pUz)) ~ 0.3. This means that the physical Nambu monopole has 
a size larger than 1 in lattice units. Therefore the two different lattice definitions locate 
it sometimes differently. 

It is clear from Fig. 5 that the percolation probability is the order parameter of the 
transition from the symmetric to the broken phase. We did not investigated the order 
of the transition. However, according to the previous investigations of the SU (2) Higgs 
model [26] we expect that for our choice of the Higgs potential it could actually be a 
crossover. 

In order to compare the position of the transition between the symmetric and broken 
phases with the point where the percolation probability vanishes, we investigated the 
susceptibility x = {H"^) — {H)"^ extracted both from Hz = J2y Z'^y and Hw = J2y |W^a;j/P- 
The dependence of x on 7 along the line of constant a = 1/128 is shown in Fig. 6. 

The magnetic energy ASp carried by a Nambu monopole is presented in Fig. 7. 
The excess of the link action near the monopole worldline ASi is shown in Fig. 8. 
The behavior of both variables show that a quantum Nambu monopole may indeed be 
considered physical object. 

5.4. Relation between the lattice model and continuum physics 

The real continuum physics should be approached along the line of constant a^, i.e., 
along the line of constant physics (at this point we omit consideration of 6w and Mh)- 
From our data it follows that the ultraviolet cutoff A = = (90 GeV)/Mz grows with 
decreasing 7 along the line of constant physics. This dependence is not far from the 
tree - level estimate A/GeV = 90/Mz = 90^/3/7 cos6'vk ~ 310/^/7. It occurs that A is 
increasing slowly along this line with decreasing 7 and achieves a value close to 430 ± 40 
GeV at the transition point between the physical Higgs phase and the symmetric phase. 
(At (3 = 15 the transition occurs at 7c = 0.92 ± 0.02, where the Z boson mass in lattice 
units is evaluated to be Mz = 0.21 ±0.02.) Therefore, we claim that up to finite volume 
artifacts the largest achievable value of the ultraviolet cutoff is around 430 ± 40 GeV if 
the potential for the Higgs field is considered in the London limit. 
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It is interesting to understand what happens with this maximal value of the 
ultraviolet cutoff if the Higgs potential would contain a finite scalar self coupling A 
[9j. Then for the case of the SU{2) x U{l)/Z2 symmetric model we chose the action of 
the form 

S = f] J2 ((l-|Trf/pCosep) + i(l-cos2^p)) + 

plaquettcs 

- 7E ^e($+f/.,e*''-<l>) + Ed"^-!' + ^d"^-!' - 1)')' (32) 

xy X 

where $ is the scalar doublet. Here 7 = 2k, where n corresponds to the constant used in 
the investigations of the SU{2) gauge Higgs model [I5l[l6lll71[l8l|l3[20l|2ll[22l|23l[2ll 
ESI [ini [m [12] ■ The weak coupling expansion in lattice theory [13] gives the prediction 
that the maximal possible ultraviolet cutoff is achieved in lattice electroweak theory at 
infinite A. Thus we expect that the largest achievable value of the ultraviolet cutoff for 
any A should be close to the value 430 ± 40 GeV calculated in our present work. One 
can compare this result with that of the previous research (see Table 1). 

We like to note, that from the point of view of perturbation theory the energy scale 
1 TeV appears in the so-called hierarchy problem [2] . Namely, the mass parameter /i^ for 
the scalar field receives a quadratically divergent contribution in one loop. Therefore, 
formally the initial mass parameter (/x^ = —X^v"^, where v is the vacuum average of 
the scalar field) should be set to infinity in such a way that the renormalized mass 
remains negative and finite. This is the content of the so-called fine tuning. It 
is commonly believed that this fine tuning is not natural [2] and, therefore, the finite 
ultraviolet cutoff A should be maintained. From the requirement that the one-loop 
contribution to is less than 10|/i|j.| one derives that A ~ 1 TeV. 

Our lattice study demonstrates the following peculiar feature of electroweak theory. 
If we are moving along the line of constant a = 1/128, then the Nambu-monopole 
density decreases with increasing 7 (for 7 > 1). Its behavior is approximated with a 
good accuracy by the simple formula: 

p ~ e2 °8-4.67. (33) 

It is worth mentioning, that the monopole density given by (133|) is expressed in lattice 
units. The density expressed in physical units may be defined as pc = where a 
is the lattice spacing. It was mentioned above that the lattice spacing decreases with 
decreasing 7. Thus the Nambu monopole density expressed in physical units increases 
with decreasing a. 

Naively one may think that the density (in lattice units) should decrease with 
increasing ultraviolet cutoff and the physical value of the density is achieved at the 
transition point. However, it occurs that the situation is inverse. Thus the density in 
physical units does not tend to a constant at a — > 0. Therefore, it is not clear how to 
extract the physical continuum density of Nambu-monopoles from the lattice data. In 
this connection it is also important to notice that we did not find any dependence of 
the monopole density on the lattice size. 
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6. On the difference between the considered lattice models 

In the previous section we have seen that there exists no difference between tlie two 
lattice models with the gauge groups SU{2) x f/(l)/Z2, and SU{2) x U{1) at realistic 
values of the coupling constants. However, such a difference clearly exists in the region 
of large coupling constants {a > 0.1), where the phase diagrams of the two models 
do not coincide. In particular, this difference can be easily seen from the behavior of 
the hypercharge-monopole density (see Fig. 9, where the dependence on j3 is shown for 
7 = 1.5). In the same figure the density of Nambu monopoles (defined through the Z' 
field) is presented as well. It is clear that the two kinds of monopoles behave differently 
in the two models. 

The explanation of this fact may be related to the possible appearance of the 
unification of fundamental interactions at the energy scale A of about 1 TeV. Namely, 
it has been shown in [36] that if TeV physics is described by a simply connected unified 
gauge group (as in Petite Unification Models [371 [3H]), then the following relation exists 
between the additional discrete symmetry of the Standard Model and the monopole 
content of the theory, which describes TeV physics: If the electroweak theory has the 
gauge group SU{2) x U{l)/Z2, then there are topologically stable monopoles in the 
unified theory, which are composed of electroweak fields (when seen from large enough 
distances). From those monopoles the hypercharge magnetic flow 



should emanate. Therefore, at low energies they may be identified with the hypercharge 
monopoles, defined in Eq. ( l3Ti) . Such objects do not appear if the gauge group of 
electroweak theory is SU{2) x f/(l). Those topologically stable objects have masses 
of the order of A/a. At realistic values of the coupling constant they could not be 
observed at low energies within the electroweak theory with the ultraviolet cutoff A, 
as their masses appear to be much larger than the cutoff. However, if one would 
imagine that the coupling constant a in our world becomes close to unity, then the 
mass of such objects becomes comparable to A. If so, their density in the case of the 
SU (2) X f/(l)/Z2-symmetric model must exceed considerably the same density calculated 
within the SU{2) x U{1) model. 

Exactly this happens in our models, where the ultraviolet cutoff A is estimated 
to be of the order of 400 GeV. The naive expression A/a for the mass of hypercharge 
monopole (in the SU{2) x f/(l)/Z2-symmetric model) gives values comparable to A in 
the region of couplings presented in the figure, where a is found to be of the order of 
0.1. It still remains larger than the cutoff, but we should remember that the classical 
evaluation of mass may be renormalized via quantum fluctuations. Thus, we can see that 
the density of hypercharge monopoles in the SU{2) x U{l)/Z2 model indeed exceeds the 
one of the SU (2) xU{l) model in the region of small (3 (the bare value of a is l/(47r/3)). 




(34) 



Namhu Monopoles in Lattice Electroweak Theory 



16 



7. Conclusions 

In this paper we investigated lattice electroweak theory numerically at realistic values 
of the fine structure constant and the Weinberg angle. We considered the potential for 
the Higgs field in the London limit, i.e., for infinite bare scalar self coupling. 

We found that the two definitions of the theory (with the gauge groups SU{2) ® 
U{l)/Z2 and SU{2) ® f/(l), respectively) do not lead to different predictions at these 
values of the couplings. However, the corresponding models behave differently at 
unphysically large values of a. The main difference is in the behavior of the so-called 
hypercharge monopoles, which would become the Z2 monopoles of the unified theory, if 
the latter has a simply connected gauge group. 

On the phase diagram of the considered lattice model a line can be drawn, where 
the renormalized fine structure constant is close to its realistic value It should be 
remembered, though, that in order to draw the true line of constant physics (where the 
Higgs mass and the renormalized Weinberg angle are constant as well) one must vary 
the bare Weinberg angle and the bare scalar self coupling. 

Our investigation of the line of constant renormalized a for the infinite bare self 
coupling of the Higgs field allows us to draw the conclusion that values of lattice spacings 
smaller than about (430 ± 40 GeV)~^ cannot be achieved in principle for this choice of 
the potential for the Higgs field (at least, for the considered lattice sizes). It would be 
important, therefore, to consider finite values of the scalar self couplings and investigate 
finite volume artifacts in order to understand whether there is a maximal value of the 
ultraviolet cutoff in the electroweak theory (that is not related to a Landau pole in 
renormalized a and A). 

Actually, we suppose that at the point of the transition the line of constant physics 
(corresponding to the Higgs phase) stops and another line of constant physics begins 
(that corresponds to the unphysical symmetric phase). Although usually the transition 
is thought to be a crossover, we see that the physical content of the theory is changed 
drastically at the transition. Namely, Nambu monopoles appear to be condensed in 
the symmetric phase. So, these objects dominate in the dynamics of this phase of 
the theory. This, of course, contradicts all observable data. Thus the points of the 
symmetric phase cannot be associated with real physics. The percolation probability 
of Nambu monopoles appears as an order parameter for this transition. Therefore we 
conclude, that the given transition may belong to the class of the transitions of the so 
- called Kertesz type (see, for example, [Ml SO]). That's why one may suppose that the 
position of this percolation transition may not correspond precisely to the position of 
the transition determined with the aid of other physical quantities. However, according 
to our results, for example, the position of the percolation transition coincides within 
statistical errors with the position of the maximum of the susceptibility extracted from 
the Higgs boson creation operator (see Figs. [5] and [6]). 

It is worth mentioning, that the gauge boson mass does not vanish at the transition 
point. So, it differs from zero in the symmetric phase of the theory, not only in 
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the physical Higgs phase. This is in accordance with the previous numerical data 
[9l[26l[T5ll6lEl[l8l[T9l[20l[2ll[23[23 For pure SU{2) theory 

(the limit 7 — 0) the situation is more complicated: The gluon propagator contains a 
mass parameter, but the dependence on the momentum is not consistent with the usual 
mass pole (see, for example, [S]). 

We have found that the Nambu monopole density on the lattice increases with 
increasing ultraviolet cutoff. Although we obtain this result for fixed bare Weinberg 
angle and infinite bare scalar self coupling, we expect that it should remain unchanged 
in the full theory with variable 6w and A. Thus we conclude, that the conventional 
electroweak theory is not able to predict the density of Nambu monopoles. This is, 
however, not a surprise because the Standard Model should be considered as a finite 
cutoff theory. According to common lore and in accordance with our numerical results 
discussed above, the Standard Model cannot describe nature at energies above 1 TeV. 
On the other hand, in [1] the Nambu monopole mass was (roughly) estimated to be 
in the TeV region. This rough estimate, however, is obtained based on the classical 
consideration of monopole configurations. In quantum theory the density of monopoles 
is defined by the balance of entropy and self energy. We found that the density of the 
monopoles approaches a value of the order of 0.1 (in lattice units) near the transition 
to the unphysical phase of the theory. Therefore, it is natural to suppose that the 
self energy of these objects decreases when the transition is approached within the 
Weinberg - Salam model. We indeed observe the decrease of the action near monopole 
trajectories (see Fig.(I7j) and Fig ([H])). So, we expect that the value of mass evaluated 
near the transition should be significantly less than the classical value calculated in [1]. 
On the other hand, it is natural to suppose that the monopole mass approaches its 
classical value deep in the Higgs phase. There the Nambu monopole currents are not 
observed in our study as it should when one looks for an object with a mass above 1 
TeV at energies of the order of 100 GeV. 

To conclude, the description of Nambu monopole physics based on the lattice 
Weinberg - Salam model with finite ultraviolet cutoff seems to us incomplete. The 
correct nonperturbative description may, therefore, be obtained within the lattice theory 
that involves the description of TeV - scale physics. Nevertheless, on the basis of 
electroweak theory only we are able to reach the conclusion that Nambu monopoles 
appear as real physical objects: Our numerical results show that the action density 
near the Nambu monopole worldlines exceeds the density averaged over the lattice in 
the physical region of the phase diagram. This confirms the classical consideration of 
[1], where corresponding classical configurations were found to carry energy. Another 
important piece of information about the structure of the lattice Weinberg - Salam model 
comes from the consideration of the percolation of monopole currents. The percolation 
probability is found to be an order parameter for the transition between the symmetric 
and broken phases. In the unphysical symmetric phase of the model Nambu monopoles 
are condensed, which means that this phase of the lattice model indeed cannot be 
associated with the real continuum physics. 
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Figure 1. The phase diagrams of the models in the (/3, 7)-plane. The dashed line 
and the dot-dashed line represent phase transitions in the SU{2) (g) U{1) model. The 
transitions for the SU{2) (g) C/(l)/Z2-symmetric model are represented by the solid 
lines. The double-dotted-dashed vertical line on the right-hand side of the diagram 
represents the line, where the renormalized a is constant and equal to 1/128. 
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Figure 2. J2x y < ^x^y > as a function of r = \xo — yo\. Here i? = 16 is the 
lattice size in "space" direction. 
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Figure 3. Potential for right-handed leptons as a function of 1/r for /? = 15, 7 = 1. 
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Figure 4. Renormalized coupling l/an as a function of f3 for 
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Figure 5. Nanibu monopole density and percolation probability as a function of 7 
along the line of constant l/aji = 128. 
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Figure 6. Susceptibility x — {H"^) ~ {H)'^ along the line of constant 1/aij 
Here Hw is the operator defined in while Hz is defined by (PO)) . 
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Figure 7. Excess of the plaquette action near monopole trajectories along the line of 
constant I/an = 128. 
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Figure 8. Excess of the link action near monopole trajectories along the line of 
constant l/ajj = 128. 
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Figure 9. Densities of hypercharge monopoles and Nambu monopoles (extracted from 
Z') as a function of /? for 7 = 1.5 for both models. 



